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We propose a chiral random matrix model which properly incorporates the flavor-number depen- 
dence of the phase transition owing to the Ua(1) anomaly term. At finite temperature, the model 
shows the second-order phase transition with mean-field critical exponents for two massless flavors, 
while in the case of three massless flavors the transition turns out to be of the first order. The 
topological susceptibility satisfies the anomalous Ua(1) Ward identity and decreases gradually with 
the temperature increased. 

I. INTRODUCTION 

Chiral symmetry breaking manifests itself in the accumulation of the Dirac modes with zero eigenvalues through 
the Banks-Casher relation[l|. In the chiral random matrix (ChRM) theory, the Dirac operator is restricted to the 
space of the constant modes and replaced with a matrix of random entities, retaining the global symmetries of QCD. 
The ChRM theory has been successful for providing a universal framework to investigate the correlation properties 
of the low-lying Dirac eigenvalues in the so-called epsilon regime in QCD@, i]. The ChRM theory can be also 
regarded as the simplest schematic model for qualitative study of the QCD-like phase diagram from the viewpoint 
of chiral symmetry Q. Effects of the external fields such as the temperature tJ^Q, the quark chemical potential [i 
0; H IE El E3| and the Polyakov loop[ll[ have been investigated within the ChRM models. 

The conventional ChRM model predicts a second-order phase transition at finite temperature T0,Q, irrespective of 
the number of flavors Nf. This is a shortcoming as a model of QCD. Based on the universality argument the chiral 
transition for Nf = 2 with the Ua(1) anomaly is expected to be of the second order with the 0(4) critical exponents, 
while the transition becomes of the first order for Nf > 3. Even in a mean-field description, the Ginzburg-Landau 
effective potential for Nf = 3 involves the Ua(1) breaking determinant term[l3|, [l4j|, which gives rise to a first-order 
transition. However, it is unknown so far how to incorporate the UA(l)-breaking determinant term into the ChRM 
models. 

The explicit Ua(1) breaking, or the anomaly, is included in the ChRM models by treating v exact zero modes with 
definite chirality, which should be interpreted as the zero modes associated with the gauge field configurations of the 
topological charge v. The fluctuation of v is effectively modeled with the instanton ensemble. The ChRM model 
supplemented with a Gaussian distribution in v has been shown to provide a screening of the topological charge 
fluctuation at zero temperature, and thus a resolution of the Ua(1) problem[l5j]. At finite temperatures, however, 
this ChRM model results in an unphysical suppression of the topological susceptibility (lH \v\ . The resolution of the 
Ua(1) problem in the ChRM model is then limited only at zero temperature. A modification of the model has been 
proposed infill \v\ to remove the unphysical suppression of the topological susceptibility at finite temperature. 

The ground state is fixed exactly with the saddle point condition in the thermodynamic limit. Because the topolog- 
ical charge v fluctuates around zero with the variance of 0(V), as dictated by the central limit theorem, the nonzero 
v configurations are irrelevant to the ground state condition in the ChRM model, and therefore the Ua(1) breaking 
seemingly cannot alter the order of the phase transition at finite temperature. 

In this paper we propose a ChRM model involving the Ua(1) breaking effect, which describes a first-order phase 
transition for Nf = 3 and at the same time reproduces the physical behavior for the topological susceptibility. The 
proposed model consists of the "near-zero" modes and the "topological zero" modes following Refs. [18|, The 
latters are interpreted as the zero modes with the right and left chiralities, respectively, associated with instantons 
and anti- instantons. We introduce number distributions of the topological zero modes based on the instanton gas 
picture. Summing over the number of the topological zero modes in the partition function, we can derive the Ua(1) 
breaking interaction. Furthermore, this model satisfies the anomalous Ua(1) Ward identity [2lj][l5[ and predicts 
physical temperature dependence for the topological susceptibility (ill [l7|. 

Interplay between the topological zero modes and the near-zero modes is studied in Ref. [2fj| to derive the instanton- 
induced interactions in the coarse-grained instanton-liquid model. The 'topological' zero modes are also introduced 
in Ref. [l8j from the determinant interaction in the Nambu-Jona-Lasinio (NJL) model in 0+1 dimensions, but the 
potential there is unbound in the Nf = 3 case. 

This paper is organized as follows. After reviewing the conventional ChRM model in the next section, we introduce 
an extended ChRM model which involves two kinds of zero modes, the "near-zero" modes and the "topological zero" 
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modes, motivated by the instanton gas picture in subsec. Ill Al Unlike in the conventional ChRM model, the total 
number of the topological zero modes is assumed to vary according to the instanton distribution. In subsec. Ill Bl 
we propose the binomial distribution for the instanton numbers in a finite space-time volume, and show that after 
summing up the topological zero modes the resulting effective potential has the Ua(1) breaking term and bounded 
from the below. General properties of the ground state and fluctuations of the model are presented for a certain set 
of parameters and equal quark masses in sec. Iim Section HVl is devoted to discussions and summary. 

II. CHIRAL RANDOM MATRIX MODEL 

The QCD partition function with Nf quark flavors of mass my is written as 

oo oo I Nf \ 

z qcd = £ e ** Z QPD = £ eWndetp + m,)) , (1) 

V— — OO U — — OO \f— 1 / y 

where (•••)[, denotes the average over the gauge field configurations of fixed topological charge is, and the 9 parameter 
has been introduced. 

The ChRM model is defined for fixed v conventionally asQ 

r *t 
Z v = d We~ N ^ wfw Y[ det(D + m f ) , (2) 
J f=i 

where the Dirac operator has been replaced with an anti-hermite matrix of constant modes 



( iW + itl N _ lu y 2 

1 iW^ + itl 



l N-\v\/2 



(3) 



with an (N + v/2) x (N — u/2) random complex matrix W in the chiral representation 75 = diag(Ijv+i,/2 7 — In—v/2)- 
Here we have introduced the effective temperature t as a deterministic part in the Dirac operator. It may be interpreted 
as the lowest Matsubara frequency t = ttT. It is readily shown that the matrix D has \v\ exact zero eigenvalues with 
definite chirality, which are interpreted as the exact zero modes associated with the topological number v. The total 
number of modes 2N is finite and proportional to the space-time volume V; N/V = 0(1). 

After bosonization of this model we obtain the effective potential fl v (S) on the chiral manifold; 

z „ . J dSe -^ Uet ((S + m/)(St +m/) + , r <4 x {*;<^);_„ < : ^ °> 

= / dSe- 2Nn ^ , (4) 



where S is an Nf x Nf complex matrix of the chiral order parameter (see next subsection). The complete partition 
function is obtained after summing over v 



•IN 



Z e = e~^^e iv6 Z v . (5) 



v=-1N 



Here we need to supplement the distribution of v characterized by the quenched topological susceptibility r, which is 
determined by the pure gluonic dynamics. Note that the range of the topological charge v is limited within ±2Af. 

In the thermodynamic limit TV — > 00, the ground state is exactly determined by the saddle point equation. Since 
the fluctuation of v scales as v 2 ~ N and nonzero v contribution to the ground state becomes negligible, the ground 
state is given by the saddle point condition for VL V (S) with v = 0, which is symmetric under U(Nf) x\J(Nf) for 
vrif — 0. Furthermore, the Nf dependence is factored out in the potential flo(S) for S oc ljy, and mf = 0,M 
yielding a second-order phase transition irrespective of Nf. Concerning the fluctuation properties, the model ([5]) 
at zero temperature gives a nonzero singlet pseudo-scalar susceptibility Xpso m the chiral limit, resolving the Ua(1) 
problem [l5|. At finite temperature, however, this leads to an unphysical suppression of the topological susceptibility 
as mentioned in Introduction. 



[*] S oc Ijvj can be derived automatically from the saddle point condition in the case of equal quark masses. 
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A. Model with near- and topological zero modes 



Let us consider a variation of the conventional ChRM model with the instanton gas picture in mind. An isolated 
instanton is a localized object accompanying a right-handed fermion zero mode. In a dilute system of N + instantons 
and N_ anti-instantons, we expect N + right-handed and N_ left-handed zero modes (even at finite temperature). In 
an effective theory at long distances, effects of the instantons should be integrated out, which will result in Ua(1)- 
breaking effective interactions 1J| . The fundamental assumption in our modeling is the classification of the constant 
modes into the "near-zero" modes and "topological zero" modes [HI]. We deal with the 2N near-zero modes appearing 
in the conventional models and include additionally the N + + A_ topological zero modes which we regard as the 
modes accompanied by the instantons. Distributions of N + and A_ lead to fluctuations of the topological charge 
v = A+ — N- as well as the total number N+ + A_ . Eventually we shall sum over 7V+ and N- with the mean value 
of O(N) and take the thermodynamic limit N — > oo. 

For definite numbers of zero modes, we write a Gaussian ChRM model as Ha 



ryN 



= / dAdBdXdYe~ N ^ AAt+BBf+xxt+YY ^l[d C t(D + m f ) 
J f=i 



(6) 



with 



D 





iA^ + itl N 




iA + itl N iX\ 

iY 

iyt o iB 

iBi / 



(7) 



Here we use the chiral representation in which 75 = diag(ljv 5 — Ijv, 1jv +5 — 1jv_)- Then the matrix D, which satisfies 
{D, 75} = 0, has a block structure with complex matrices A, B, X and Y as shown above. The N x N complex matrix 
A acts on the near-zero modes, while the rectangular complex matrix B of size N + x A_ acts on the topological zero 
modes. The matrices X € £NxN- an( j y g ^NxN + re p resen t the interactions between the near-zero modes and the 
topological zero modes. In the following we take the quark mass mj of flavor / to be diagonal cx l2./v+./v + +./v_ in the 
space of the zero modes. Notice that the temperature term t is introduced only for the near-zero modes, while the 
topological zero modes are assumed to be insensitive to the temperature t. This discrimination is physically legitimate 
in the instanton gas picture where each topological zero mode is localized around an (anti-)instanton and its eigenvalue 
is not much affected by the anti-periodic boundary condition in the temporal direction unless the temporal size 1/T 
becomes the same order of the typical temporal extent of the instanton or the topological zero mode. I**' 

The effective potential on the chiral manifold can be derived from eq. ([5]) in a standard manner. First, we recast 
the partition function in the form of the integration over Grassmann variables of the near-zero modes ip and of the 
topological zero modes x : 



dct(D + rrij 



dtp^ dipdx d\ 



exp 



(4,Vi,xLxii)(£> + m/) 




(8) 



where, for each flavor, ip^ L and i^r.l have N components while x1(_r) an d Xr(l) have components. Then 

integration over the random matrices A, B, X and Y with the Gaussian distributions gives rise to the four-fermion 
interaction term. We unfold this term with the aid of the Hubbard-Stratonovich transformation as 



exp 



I 



NY? 



dS exp 



xl f X R )Sfg 



xUx 9 L )S) g ] exp [-A£ 2 tr£tS] 



(9) 



with an Nf x Nf complex matrix S. Here summation over the flavor indices / and g = 1, • • • , Nf should be under- 
stood, and the zero-mode sum is also implicit in (e.g.) %l>lfip 9 R = Y^^u^Ri' This transformation allows us to perform 



[**] If T is high enough to modify the nature of the topological zero modes, one should deal with the modification of the instanton 
configurations and distributions, too, which is beyond the scope of this study. 
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the Grassmann integration. We thus obtain the desired form for the partition function: 

Z$ +tN _ = [ dSe~ NsHTSts det N [(S + M)(S^ +M^)+t 2 ] det N +(S + M) det N -(S^+M^). (10) 
The mass matrix M. = diag(mi, • • • , mjv,) will be taken as a general Nf x Nf matrix of the source field later in 

sec. inn 

In the chiral limit M. = 0, the integrand of (p~0|) is invariant under S\J(Nf) x SV(Nf) transformation, S — > USV" 1 
with U, V £ SU(Nf), in addition to Uy(l)- On the other hand, under Ua(1) transformation, S — > e S, it acquires 
a phase factor e 2lAr / Sl/ due to the difference v = N + — iV_ in the powers of the determinants. In other words, the 
nonzero v breaks the Ua(1) symmetry explicitly down to liWj- 



B. Distributions of topological zero modes 

The complete partition function is obtained after summing Z$ n_ over the 'instanton' numbers N + and AL_ 
with an appropriate weight reflecting the pure gluon dynamics.^' Although the distributions of N + and iV_ will be 
correlated in general, here we assume independent distributions P{N±) for N+ and N-, i.e., 

Z e = eW6p ( N +) p ( N -) Z N + ,N_= / dSe- 2Jvn(S;t - m ' e > . (11) 

N+,N- 

First we consider P(N±) in a dilute instanton gas picture. For one-instanton configuration, one may assign a weight 
k compared with a no-instanton configuration, and multiply a factor N oc V taking into account the integration over 
the instanton location. For a configuration with N + instantons, we have a weight 

Ppo(N + ) = ^(*N) N + (12) 

where N + l is the symmetry factor [TJ]. The same distribution is assumed for iV_ as well. This is nothing but the 
Poisson distribution up to a normalization. The summation with Pp Q (N±) in eq. (|11[) results in the exponentiation 
of the determinant term(l8l|: 

J2 Po =iE 2 trS'S' t - i lndct [(S + M)(S* + M^) + t 2 ] - \n\e w det(5 + M) + e~ ie det(S' t + M^)]. (13) 

The same determinant term as in flp Q is commonly incorporated in other effective models as the "anomaly term" to 
break the Ua(1) symmetry. M In the ChRM model, however, this potential is unbound for Nf = 3. Indeed, for large 
S = 4>lN f the term det(S' + M.) ~ cj> 3 dominates over the other terms in f2p . 

It should be noticed here that the fermion coupling distorts the N± distribution itself. With including the deter- 
minant term of the topological zero modes in eq. (fTDj) . the effective distribution for N + reads 

P Po (N + ) = {ltNd f+ (14) 

with d = \det(S + M)\, and similarly for 7V_. We find that the average value of N± increases indefinitely with 
increasing d as (N±) = nNd. However, the possibility of infinitely many topological zero modes N± should be 
avoided in the ChRM model as a low-energy effective theory within a finite volume V oc N. 

Here we regularize the distribution by setting explicitly a maximum value of O(N) for N±. We split the finite 
space-time volume into jN cells with 7 being a constant of 0(1), and assign a probability p for a cell to be occupied 
by a single '(anti-) instanton' and (1 — p) for a cell unoccupied. We exclude the possibility of the double occupation 
of a cell, which is justified by a repulsion between the instantons. Just like in the lattice gas model in statistical 
mechanics, this assumption results in the binomial distributions for N±: 

nw±) = (3£) p n hi-pv n - n± - (15) 



[f] We keep here the the number of near-zero modes 2N fixed. In Ref. |19| , the total number of the zero modes, 27V + 7V+ + was fixed. 
[J] Although one assumes no mass term in the determinant interaction conventionally, there is an ambiguity [23] ■ Our choice of the mass 
dependence has been fixed in eq. SEl . 
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For a small p and a large 7 AT, the binomial distribution P(N±) is accurately approximated with the Poisson distribution 
with the mean ^Np. But it cannot for a large p. The binomial distribution provides a stringent upper bound 7-/V 
for the number of topological zero modes N±, in contrast to the Poisson distribution. The corresponding effective 
potential for S is found to be 

fl(S;t,m,e) =|S 2 tr55 t - |lndet [{S + M){S^ + M ] ) + t 2 ] 

- ±7 [In (e ie a det(S* + M) + l) + In (e~ l£ 'a det(5 f + M r ) + l)] (16) 

with a = p/(l — p). This is the ChRM model that we propose and analyze in this paper. 

The effective potential (|16|) is bounded from the below by the ti SS> term in contrast to flp in eq. (fl"3|) . The anomaly 
terms are accommodated under the logarithms in the square bracket. Moreover, for small value of |adct(S' + m)|, it 
reduces to the potential (fl"3|) with k — -fa. On the other hand, the Poisson approximation fails when (N±) becomes 
0(7 AT). We stress here again the fact that the distributions P(N±) is deformed in the presence of the coupling with 
the topological zero modes as 

p(n + ) = ( ( P d) N +(i- P y N - N +, (17) 

which changes the probability p to an effective one p = pd/{pd + 1 — p). Similarly for A^_. Accordingly the mean 
number of the zero modes is modified to (N±) = jNp. 

Two remarks are here in order. In the conventional models, the total number of the modes is kept fixed, while the 
fluctuation of v is allowed for resolving the Ua(1) problem. On the other hand, in the instanton gas picture, both 
of N + and AT_ are expected to fluctuate naturally. As stressed in this subsection, the mean number of the modes 
depends on the magnitude of S, which gives an O(N) effect. The anomaly term appearing in the effective potential 
in turn affects the saddle point value of S. This point is essential to the first-order transition for Nf = 3, which is 
overlooked in the earlier works. 

Secondly, in cq. (|16[) one needs to introduce a dimensionful scale to compensate the dimension of determinant if the 
dimension of mass is assigned to M. as well as S. Here we leave all the quantities dimensionless for demonstration of 
general feature of the model. 



C. 9 dependence 

The variance of the topological charge v = N+—N- for the binomial distribution is computed as 2Nt = 2N~fp(l—p), 
where r is the quenched topological susceptibility. In the presence of the fcrmion coupling, this susceptibility will be 
replaced with 



7P(1 -P)=J-, 



ad 



(ad+l) 2 ■ 

We can confirm this fact also by rewriting the anomaly term in eq. (|16|) as 



1 + \a dct(S* + M)\ 2 + 2|adet(5 + M)\ cos id - - In 



det{S + M) 
det(£t +JWt) 



(18) 



(19) 



from which we find again the replaced topological susceptibility f as a coefficient of 9 2 . Such a series in of the 
anomaly term gives a connection to the general arguments on the 77' in the 1 /N c expansion [23[ . 

To make a clear connection with the conventional model, let us ignore for a moment the fluctuation of the total 
number N+ + N_ and apply the Gaussian approximation for the v distribution. We then obtain the partition function 

as 



=const. 



dS e 



-ATS 2 trS t S 



det^ [(S + M)^ +M f ) + t 2 



x det^ +N -^ 2 [(S + M)(&+Ml)] 



exp 



—Nt 



_i_ det(5 + M) 
2 n det(5t + A<ft) 



(20) 
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This is almost identical with the model discussed in Ref . Jig, [13] , which reproduces the screening of the (unquenched) 
topological susceptibility (see eq. as measured on a lattice at finite temperatures. However, this model ([2D)) fails 

to describe a first-order phase transition at finite temperature for Nf=3. In fact, the anomaly term appears only as 
a phase in (|20p in contrast to (|19[) and drops out when we determine the magnitude of S = S* in the ground state 
(with M = A4'). The variation of N + + is essential for the anomaly term to affect the saddle point condition, 
and thus the order of the phase transition. 

III. GROUND STATE AND FLUCTUATIONS 

In this section we shall study the ground state properties of the system with equal mass, Ai = mljv,, for simplicity. 
Setting S = 4>\n s with real </> and with 9 = 0, we obtain a simple form of the grand potential: 

Q((f>;t,m) =\N S Y?4> 2 - iiV/ln [(0 + rn) 2 +t 2 ] - 7 m |a(0 + ro)^ +l| . (21) 
The factor Nf cannot be factored out in the potential f2 because of the anomaly term here. 

A. Chiral phase transition: ground state 

In the thermodynamic limit, N — * 00, the ground state can be analyzed with the solution of the saddle point 
equation 

2 _ + to _ a^ + m^f- 1 _ 
9 {^ + m) 2 +t 2 1 a(<f> + m) N f + 1 [ ' 

The scalar quark condensate is related to the solution as 

(M) = 1 — ^- In Z(m) = -£ 2 0. (23) 

2NN f dm y 1 v 17 



Without the anomaly term a = 0, eq. (|22|) recovers the flavor-independent gap equation of the conventional model, 
which has the so 
leading order as 



which has the solution 0q = S 2 — t 2 in the chiral limit. One can estimate the effect of small m and a on tp 2 at the 



At lower temperatures 2t 2 < S -2 the value of the condensate is decreased by the quark mass term, on the contrary to 
our intuition. This is because the leading order term of m appears in a combination — 2to</>/ {(f) 2 + 1 2 ) in the potential, 
which simplifies to — 2m/ (f> for t = 0, favoring smaller <fi for the potential to be more stabilized. 

There is no chiral symmetry for Nf = 1 because of the anomaly term. Let us discuss phenomenologically interesting 
cases, Nf = 2 and 3. 

1. N f = 2 

Expanding the potential Q around = 0, we find 

n = c + c 2 <j) 2 + c 4 <?!) 4 -h(j) + 0{<f, 3 m, to 2 ) (25) 

with Co = — lnt 2 , C2 = S 2 — — t~ 2 , C4 = (7a 2 + 1~ 4 )/2 and h = 2(aj + t~ 2 )m. This is the standard form of the 
Landau-Ginzburg potential for a second order phase transition (to = 0) with the critical temperature t c 

(26) 



v/s 2 - 07 

Inclusion of the anomaly term increases the value of t c . The behavior in the vicinity of the critical point is characterized 
by the mean-field exponents; we find (3 = 1/2 and <5 = 3, respectively, from the solutions cj> 2 = et~ 2 /c4 cx e 2 ^ for 
e = (t c - t)/t c with < e < 1, and <f> = [toE 2 /(2c 4 )] 1 / 3 cx m 1 ^ with to ^ at t = t c . 
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2. N f = 3 

Expanding the potential Q around (f> = 0, we find 

n = c + c 2 (f> 2 + c 3 <p 3 + c 4 4 - H + 0(4> 6 J 2 m, to 2 ) (27) 

with c = -(3/2) lnt 2 , c 2 = (3/2)(£ 2 - i~ 2 ), c 3 = -cry, c 4 = (3/4)t~ 4 and h = M~ 2 m. The U A (1) symmetry is 
explicitly broken even for to = due to the anomalous 3 term, which leads to a first-order phase transition. 

In Fig.[TJ we display the chiral condensate ^ as a function of the temperature t and the quark mass to. For numerical 
demonstration, we have chosen the parameters as £ = 1, 7 = 2 and a = 0.3 in this paper. We clearly see the second 
order transition for Nf = 2, while the first order transition for Nf = 3 in the chiral limit to = 0. In Nf = 3 case, as 
we increase the current quark mass to, we find a terminating point of the first order line at m c = 0.0265. 

For Nf > 4, the anomaly term only affects the coefficients of cf> n (n > Nf) in the series expansion of the potential. 



B. Mesonic masses 



Once the ground state acquires the nonzero expectation value S = 4>1nj 7^ 0, the symmetry is spontaneously 
broken down to XJy(Nf); S — > USU~ 1 . Degeneracy of the vacua in the chiral limit dictates the massless fluctuations 
S = 0e l7r A /(^"W ~ (f) -f in a X a /y/2 corresponding to the generators of the broken symmetry. 

To find the mesonic masses, we use the parametrization S = </> + A a (er a + i7r a )/v / 2 with <r a , ir a G R and the V(N) 
generators A a normalized as tr[A a A fc ] = 25 ah (a = 0, • ■ • , iV 2 — 1). Using a formula for a matrix X with a small 
parameter e, det(l + eX) = 1 + etrX + ^e 2 [(trX) 2 — trX 2 ] + C(e 3 ), we can easily expand f2 ([16]) around the saddle 
point solution to define the massses with Cl = fl Q + ^M 2 a u a2 + ^M 2 sa ir a2 + ■ ■ ■ . 

The flavor non-singlet masses in the scalar and pseudo-scalar channels, respectively, are found to be 

M 2 =T 2 - 1 . 2^ + to) 2 a (0 + TO)^- 2 

((/) + m) 2 +t 2 [{(/) + m) 2 +t 2 } 2 l a((j> + m) N f + 1 1 1 ' 



1 a((j) + m) N f- 2 

< + m) 2 + t 2 1 a((j) + m) N f 



»r2 V 2 * -r ^ 

M ps = S - n , I I2 - 7 „,^ , _, N , , -, • (29) 



With the saddle point equation, we can re-express the it mass as 

M 2 ps = (30) 

y <fi + TO 

In the chiral limit, we have the massless it. With the explicit breaking of to, we find a relation reminiscent of the 
G ell-M ann- O akes- Renner , 



((f) + m) 2 M 2 s = mS 2 (0 + to) ~ —m^ip), 



(31) 
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FIG. 2: Temperature dependence of the mesonic masses in the flavor-singlet scalar (sO) and pseudo-scalar (psO) channels and 
in the flavor-nonsinglet scalar (s) and pseudo-scalar (ps) channels in the chiral limit (m = 0; thick lines) and with explicit 
breaking (m 7^ 0; thin lines) for Nf — 2 (left) and Nf = 3 (right). As nonzero quark mass, we set m = 0.1 for Nf=2 and 
m = m c = 0.0265 for Nf = 3. 



if we identify + m as the pion decay constant f^. Interestingly, there is a mass hierarchy, M 2 : E 2 : M 2 S = 2:1:0, 
in the chiral limit at zero temperature t = 0. 

On the other hand, the flavor singlet masses for the scalar and pseudo-scalar singlet channels have an additional 
contribution from the anomaly term as 

Ml =M 2 - AMq , (32) 
Af 2 s0 =Ml a + AM 2 (33) 



with 



AM = Nfj- 



[a((j) + m) N f + 1] 



r 



= N fTl so ■ ( 34 ) 

The would-be Nambu-Goldstone (NG) mode in the singlet channel becomes massive due to the coupling to the Ua(1) 
interaction AMq . This same effect appears as a reduction in the scalar singlet channel. This mass gap is related to 
the (replaced) quenched susceptibility f, similarly to the Witten-Veneziano formula[23j|. It is interesting to note that 
for Nf > 3 the Ua(1) breaking effect disappears AMq = from the meson masses when the condensate vanishes in 
the symmetric phase in the chiral limit [24l] . In contrast, iV/-point functions are affected by the Ua(1) breaking term. 

In Fig. [5] we show the mesonic masses in the flavor-singlet scalar and pseudo-scalar channels and in the flavor- 
nonsinglet scalar and pseudo-scalar channels. Both for Nf=2 and 3, the flavor-singlet pseudo-scalar meson acquires 
the nonzero mass via (|33[) . 

C. Susceptibilities 

The scalar and pseudo-scalar susceptibilities are defined as the responses to the external fields M. = (s a + ip a )A a /v / 2 
in eq. (fTtj|l: 



B 2 



2 



(35) 
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In our basis, those susceptibilities are found to be diagonal xi \ s ^ ^ ab with a simple form 



X =X 



(0). 



1-X (0) / S2 



,(o). 



M 2 



(36) 



with = S 2 — M 2 for each channel. A susceptibility diverges when the corresponding M 2 vanishes (e.g.) at a 
critical point where the effective potential £l(S) has a flat direction. Especially, in the broken phase, the susceptibilities 
of the non-singlet NG modes are indefinite. 

The topological susceptibility, x t op is defined as the response to the angle 9 



Xtop 



2N d9 2 



InZfi 



(37) 



)=0 



Using the saddle point solution with small 9, parametrized as S(9) = (f> + i?7o(#)A /\/2, the susceptibility in the mean 
field approximation is obtained by differentiating fl(S(9); 9) around 9 = as 

_ d 2 n / 8 2 ii \ 2 (d 2 n 



_^ N f t 2 1 



+ mf M 2 



psO 



mY, 2 (<j) + m) 



(38) 



Xtop 



(39) 



Here f is the modified susceptibility defined in cq. (|18[) in the previous section. As is well known|25|. the most 
prominent effect of the fermion coupling is the screening of the topological susceptibility via the contribution 



T, 2 m((f> + m) 
N~f 



Ml s {4> + m.f 



(40) 



In the quenched limit, Nf — > 0, the x t op recovers the quenched susceptibility r = 7p(l — p), while in the massless 
limit the susceptibility Xtop is screened completely to zero. 

By changing the variable S — > S = Se ie / Ni under the integral $W\i and then applying the saddle point approxi- 
mation, we find £l(t,m,9) = Q(t, me ie / Nf , 0). Here the 9 appears only in the combination with the mass m. In this 
latter form, the source term becomes a linear combination of the scalar and pseudo-scalar sources in the flavor-singlet 
channel: M = me i0 / N f = (s + ip )\ Q /y/2, and the Ua(1) relation for the topological susceptibility [2l| is immediately 
derived as 



Xtop — r) AT 




dlnZg d 2 s 



)=0 



(41) 



We have seen that the pseudo-scalar meson in the flavor singlet channel has nonzero mass ([3"3"| because of the Ua(1) 
breaking term, and accordingly the pseudo-scalar singlet susceptibility remains finite in the broken phase in the chiral 
limit. Thus for the small but nonzero quark mass m, the topological susceptibility Xtop decreases following the chiral 
condensate (tpip) ~ <p with increasing temperature t. 



IV. DISCUSSIONS AND SUMMARY 



In this work we have considered the ChRM model with the near-zero and topological zero modes. It was known 
before that summation over the topological zero modes with the Poisson distribution results in the determinant 
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FIG. 3: Temperature dependence of topological susceptibilities for Nf=2 (left) and 3 (right). 



interaction, but unfortunately it gives a pathologic unbound potential for the ChRM model with Nf = 3 IS]. On 
the other hand, the Gaussian distribution for topological charge v with the total number of the modes fixed, leads 
to the log-determinant type interaction (1201) . This term in the potential resolves the Ua(1) problem at zero and 
finite temperatures as suggested in Ref. [161 Il7l ] , but still yields a second-order phase transition at finite temperature, 
irrespective of Nf. 

We have proposed that (i) the numbers of the topological zero modes with right and left chiralities, respectively, 
have the distributions related to instanton dynamics in a finite volume, and thus (ii) these distributions have an 
upper bound of O(N). We have adopted the binomial distribution as such a distribution. This gives rise to the stable 
potential which describes the chiral phase transitions of the second and the first order depending on the number of 
flavors Nf=2 and 3, respectively, and resolves the Ua(1) problem as well. We have confirmed through numerical 
evaluations that the proposed model also reasonably reproduces temperature dependence of meson masses, (pseudo- 
) scalar susceptibilities and topological susceptibility. Notably, the topological susceptibility is consistent with the 
universal quark mass dependence, satisfying the anomalous Ward identity. 

In the conventional ChRM models, the near-zero modes themselves are customarily assumed to emerge from the 
instanton dynamics, in contrast to the explicit separation of the near-zero modes and the topological zero modes in 
our modeling. To assess the foundations and relations of these models, a rigorous analysis of these low-lying modes 
based on the microscopic dynamics would be needed, which is beyond the scope of our current study. 

We have adopted the independent binomial distribution for N± as a simple improvement from the Poisson distri- 
bution in order to obtain a stable effective potential. However, the distributions of N± can be correlated non-trivially 
in general. It will be interesting to examine other possibilities for number distributions of the topological zero modes. 
Incidentally, we note here that the determinant interactions of the from Q13p in the NJL models in 3+1 dimensions also 
lead to the unbound effective potentials for a large quark condensate beyond the cutoff although the local minimum 
is usually chosen as the physical ground state. 

Concerning phenomenological applications, first we need to tune the model parameters S, a and 7 as well as the 
light and strange quark masses m uc i and m s so as to reproduce the empirical properties in the vacuum. Extension 
to the case at finite baryo-chemical potential is straightforward [J H|, which allows us to study the phase diagram of 
the ChRM model with the flavor dependence in the space of the temperature, the chemical potential and the quark 
masses. Especially the existence of the analog of the QCD critical point (s)[26j will be an important subject to be 
studied. Application of this model at finite isospin and strangeness chemical potentials are also planned. Progress in 
this direction will be reported elsewhere. 
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APPENDIX A: TRINOMIAL DISTRIBUTION 

Here we deal with the trinomial distribution for N + and iV_ , as a simplest example of the non-factorizable distri- 
bution: 

p ^ - m £\-»-v ^~-e-*+-*-r~-~-" - ■ (ad 

where is the probability for a single cell to be occupied by an (anti-)instanton. Note that the distribution 

is symmetric under the exchange of + with — . Replacing P(N + )P(N_) in eq. (|lip with P(N+, N_) and setting 
p + = p- = p, wc obtain the effective potential 

n Tri (S; m, 6, t) = \Y?trS ] S - § In dot [{S + M)(S^ + M ] ) + t 2 ] - ±7 In [ae lf> det(S + M) + ae~ ie dct(S' t + M 1 ) + l] 

(A2) 

with a = p/{l — 2p). This is quite similar to eq. (|16[) . and qualitative features of the model are unchanged. The 
quenched topological susceptibility, or the variance of v = N+ — N- , for eq. (|Alj) is computed as 

t = 7p . (A3) 
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